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Abstract 



In this paper we revisit the problem of finding hidden symmetries in quantum mechanical 
systems. Our interest in this problem was renewed by nontrivial degeneracies of a simple spin 
Hamiltonian used to model spin relaxation in alkali-metal vapors. We consider this spin Hamil- 
tonian in detail and use this example to outline a general approach to finding symmetries when 
eigenvalues and eigenstates of the Hamiltonian are known. We extract all nontrivial symmetries 
responsible for the degeneracy and show that the symmetry group of the Hamiltonian is SU{2). 
The symmetry operators have a simple meaning which becomes transparent in the limit of large 
fSj ' spin. As an additional example we apply the method to the Hydrogen atom. 
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g ; 1 Introduction 

The close connection between symmetry and degeneracy has been explored since the very foundation 
of Quantum Mechanics. Famous examples include degeneracies of the energy spectra in angular 
momentum in the 3d harmonic oscillator [1] and in the Hydrogen atom (well known as the accidental 
degeneracy) [2]. Here we undertake a detailed study of the connection between symmetries and 
degeneracies of a Hamiltonian that describes the exchange interaction of two spins, S and K, and also 

3: 

H{x) = x(K + 1/2)S Z + K- S (1) 



includes Zeeman splitting for spin S: 



where S = 1 and K is arbitrary. 

Our interest in this system was motivated by experiments [3] on spin relaxation in polarized alkali- 
metal vapors. In this case S has a meaning of the total electronic spin, K is the nuclear spin, and 
the dimensionless constant x represents the magnetic field. A diagram of energy levels of Hamiltonian 
(1) for a typical value of K = 2 is shown on Fig. 1. In addition to degeneracies at x = (no 
Zeeman splitting) the spectrum displays less trivial (2K + l)-fold degeneracies at x = ±1. These 
degeneracies show up as resonances in the spin relaxation rate and provide the key evidence for a 
particular mechanism of spin relaxation (see [3] for details). The degeneracies at x — ±1 have been 
also discussed in [4]. 

Despite the simplicity of Hamiltonian (1), it is not trivial to determine its symmetry at x = 1. To 
the best of our knowledge there is no general textbook algorithm for extracting symmetries. In this 
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paper we propose such an algorithm. We demonstrate how, using only a minimal intuition into the 
reasons for degeneracy, one can find hidden symmetries whenever the eigenstates and the eigenvalues 
of the Hamiltonian are known. 
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Figure 1: Eigenvalues of H(x) = x(K + 1/2)S Z + K • S for S = 1 and K = 2. Note the level crossings 
at x — 1. 

Before proceeding further let us formulate the general problem of identifying the symmetry re- 
sponsible for a particular spectral degeneracy. Symmetries manifest themselves through operators 
that commute with the Hamiltonian, H . Consider a set of these operators A = {Aj}. If not all of the 
operators Aj commute with each other, this symmetry implies some spectral degeneracy. Indeed, if 

[A i ,H] = [A j ,H]=0 and [A^A^O, (2) 

there exists an eigenstate of the Hamiltonian \a) for which A\\a) ^ Aj\a). This state is necessarily 
degenerate, since according to equation (2) A\\a) and Aj\a) are eigenstates of H of the same energy 
as \a) (see e.g. [1]). 
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Is a particular set of operators A sufficient to account for given spectral degeneracies? The positive 
answer to this question implies that 

1. all degenerate states can be labelled by eigenvalues of a number of mutually commuting operators 
from the set A 

2. any degenerate state can be obtained from any other state in the same eigenspace by a repeated 
action of operators Aj or their linear combinations. 

Provided that conditions 1) and 2) are met one can complete the analysis by clarifying the physical 
significance of operators {A{\. An additional, more abstract, question one might ask is what is the 
group generated by operators {Ai}. 

Consider for example the degeneracy of Hamiltonian (1) at x — 0. One can introduce G = K + S 
and identify, e.g., A\ = G x and A 2 = G y . The commutator of A\ with A 2 yields A 3 = G z . Eigenvalues 
of G z and G 2 can be used to label the degenerate states, while linear combinations of G x and G y , 
G + and GL, connect all degenerate states in a given eigenspace. Finally, any operator that commutes 
with H(0) can be written in terms of G x , G y , G z , and G 2 . The group generated by the components 
of G is SU (2) and the group elements have a meaning of rotations in 3d space. 

Frequently, we encounter a situation when the spectrum is known exactly and yet the symmetries 
are hidden. In such cases the following formal expression for commuting operators serves as a useful 
starting point. One can show quite generally that, up to an operator that annihilates all degenerate 
states, all integrals of motion must be of the form 

X = Y J PnYP n (3) 

n 

Here P n is an operator that projects out the nth degenerate eigenspace and Y is an arbitrary operator. 
We prove equation (3) in the Appendix. The problem of finding the symmetry thus reduces to making 
appropriate choices for Y. This choice can be either guided by limiting cases, e.g., the classical limit, 
where the symmetry is simple to identify, or by an intuition as to what type of conservation laws 
(e.g. scalars or vectors or etc.) one expects to find. However, examples we considered show that 
many apparently different choices for the operator Y produce equivalent conservation laws. Thus, it 
is usually sufficient to explore the simplest options - the basic operators of the problem. For example, 
if the Hamiltonian is written in terms of f and p, natural choices for Y would be f 2 and p 2 if one is 
looking for a conserved scalar or f and p if a conserved vector is expected. 

The paper is organized as follows. First, we determine and analyze in detail the spectrum of 
Hamiltonian (1) at x — 1 (Section 2). We show that, apart from the zero field degeneracy at x — 0, 
there are degeneracies only at x = ±1, the spectrum at x = —1 being simply related to that at 
x — 1. Given the eigenstates, integrals of motion (3) can be evaluated without knowing projectors 
P n explicitly. In this case the general strategy is to compute the matrix elements of X in equation 
(3) and use them to determine X in terms of basic operators of the problem. We will illustrate this 
approach in Section 6. 
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In Section 3 we use (3) to explicitly derive commuting operators that connect all degenerate states 
at x — 1. We demonstrate that the symmetry group of Hamiltonian (1) at a; = 1 is 577(2). The 
physical meaning of symmetry operators is clarified in Section 4, where we consider the limit of large 
K. 

In Section 5 we rewrite the Hamiltonian in the form that makes the analogy to the large K limit 
particularly clear and allows us to establish some additional properties of the energy spectrum. In 
Section 6 we provide another example of the general approach of finding symmetries based on (3) by 
deriving the Runge-Lenz vector from (3). 



2 Energy Spectrum 



Here we analyze the spectrum of Hamiltonian (1) at arbitrary x and consider in detail the spectrum 
at x — 1. The results discussed in this section were originally derived in [6]. We will follow [6] closely 
adopting in most cases the same notation. 

First, we note that the z component of the total spin is conserved. 



[G z ,H(x)] = 



(4) 



Since S = 1, there are at most three independent states for each eigenvalue m of G z . In the basis 
\K Z , S z ) these states are \m — 1, 1), \m, 0), and \m + 1, —1). The block of the Hamiltonian (1) that 
corresponds to m for |m| < K — 1 is 



y/2 



V2(x(K + l/2)+m-l) yJ(K + m)(K-m + l) 
'(K + m)(K-m+ 1) 

y[{K 







j{K-m){K + m + \) 
m)(K + m + l) V2(-x(K + 1/2) - m - 1) 

with the eigenvalue equation 

E 3 + 2E 2 - a(x)E - K 2 - b(x) = 



(5) 



a( x ) = (1 + x 2 )(K + 1/2) 2 + 2x(K + l/2)m - 5/4 (6) 

b(x) = K 2 + K + x(K + l/2)m 

In addition to 3 x 3 blocks (5) the Hamiltonian also has two 2x2 blocks for |m| = K and two lxl 
blocks for \m\ — K + 1. Energies for these values of m are also solutions of (6). 

Degeneracies occur only at x = and x — ±1 (see Fig. 1). The energy spectrum at x = —1 is 
identical to that at x — 1 while the eigenstates are related via a unitary transformation (a rotation by 
7i around any axis in xy plane). Therefore, it is sufficient to consider only x — 1. 

The spectrum at x — 1 has the following features: 
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• There are 2K + 1 degenerate states \Dm) ("D" for degenerate states) with the energies 
Eom = —1/2 

• There is a gap at m = —K in the values of m that can be assigned to degenerate states \Dm), 
m = K, K - 1, . . . , -K + 1 and -K - 1 

• There are 2K + 2 non-degenerate states \Tm) ("T" for top states) with energies Erm > — 1/2 
and m = K + 1, K, . . . , -K 

• There are 2K non-degenerate states \Bm) ("B" for bottom states) with energies Esm < —1/2 
and m — K — 1, K, . . . , — K 



• Energies of the top \Tm) and bottom \Bm) states are related by E Tm + E Bm = —3/2 (|m| < K) 

Finally, we can use (5) to compute the energies at x — 1 and the wave-functions \Dm) of the 
degenerate states 



E Dm = -l/2 E Tm = -3/4 + ^/B(m) 2 + 1/16 £ Bm = -3/2 - # 



Tm 



(m — 1, l|Dm) 
(m, 0\Dm) 
(m + 1, — l\Dm) 



where 



V2B(m] 

B(m) = (2K + 1)(AT + m + 1/2) 



- x / (K + 1 + m)(K + 1 - m) 
2(K + l + m)(X + m) 
^/(iir + m)(A:-m) 



(7) 
(8) 

(9) 



3 Symmetries 

In this section we derive integrals of motion responsible for degeneracies at x — 1. We write them in 
terms of components of S and K and show that the symmetry group is SU(2). 

Our starting point is the general expression for commuting operators (3). Hamiltonian (1) has only 
one degenerate subspace at x — 1, so there is only one projection operator in (3). 

x = PyP (io) 

Even though we do not need to know projection operators explicitly to evaluate (3), in the case of 
Hamiltonian (1) the projector P has a simple meaning in the large K limit, so we derive it bellow 
from the eigenvalue equation (6). 

When the energy E is replaced by H and m is replaced by G z , eigenvalue equation (6) translates 
into a cubic identity for the Hamiltonian [6]. At x — 1 this identity can be factored out as follows 

(H + 1/2) \B(G Z ) 2 -(H+ 1/2)(H + 1)1 = (11) 
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where 

B(G Z ) = (2K + l)(K + G z + l/2) (12) 

States with E — —1/2 are degenerate with respect to the eigenvalues of G z . The operator in the 
square brackets, IT, is then an "unnormalized projection operator" that projects out the degenerate 
subspace. 

fl\Nm) = U\Dm) = B(m) 2 \Dm) (13) 
The operator II can be normalized to a usual projection operator for the degenerate subspace [6]. 

n (h + i/2)(h + i) , x 

= 1 - - — J A — (14) 



B(G z y B(G z y 

Now we have to identify appropriate choices for the operator Y in equation (10). Let us first try 
the simplest options - operators that are linear in components of S and K. Consider, for example, a 
conserved vector 1 

L = PKP (15) 

Let us determine the action of components of L on the eigenstates of Hamiltonian. Since P is a 
projector for the degenerate subspace 2 , 

Li\Nm) = i = 1,2,3. 

To calculate the action of Lj on degenerate states recall that due to selection rules (see e.g. [1]) the 
components of a vector can have nonzero matrix elements only for transitions m — > m and m — > mil. 
We note also the following relation between matrix elements of operators X and Y in equation (10) 

(m'\X\m) = (Dm'\Y\Dm) (16) 

where \m) and \m!) are any two eigenstates of G z with eigenvalues m and w! respectively. Using 
equations (16) and (8), we obtain 



L + \Dm) 



(K*-m>)((2K + 2m + 3) 2 - 1) 
\ {2K + 2m + 2f-l \ V > m + i ) t U > 



The action of L_ can be determined directly from equation (17) and the action of L z from equations 
(16) and (8) with m! = m. 

Note that equation (17) is not well defined for m = —K — 1. For the state \D, —K — 1) instead of 
equation (17) we get 

L±\D, -K - 1) = (18) 

1 Projection operator P is a scalar, because the Hamiltonian is a scalar. To see this the operator S z in Hamiltonian 
(1) should be written as S • b, where b is a unit vector along the magnetic field. By the same argument G z = G • b is 
a conserved scalar. Therefore, L is a vector. 

2 We use the usual notation Li = L x = L+ ~^ L ~ , L 2 = L y = L+ ~ L ~ , and L 3 = L z . 
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Equation (18) reflects the existence of a gap at m = —K in the eigenvalues of G z assigned to degenerate 
states. An operator that connects the state \D, —K — l) to the rest of degenerate subspace must change 
m by at least two. According to selection rules such an operator is neither a scalar nor a component 
of a vector. Let us first consider values of m ^ —K — 1 and revisit the problem of the gap at the end 
of the section. 

Since equation (17) is not of the standard form for a raising operator in su(2), commutation 
relations for L± are "deformed" versions of su(2) commutation relations. However, one can derive 
usual su(2) operators from L± and G z . Define 

N + = (AL)t = A(G Z )A(G Z + 1)L+ N z = (G z - l/2)P (19) 

where 

B(G« + l/2) 

Operators iVj form an su(2) algebra: 

[N i ,N j ]=e ijk N k 

All degenerate states with m ^ —K — l are connected by operators N±. Each of these states is labelled 
by N z and by an eigenvalue of the Casimir operator 

3 

N 2 = N- = N-N+ + N 2 Z + N z 
i=i 

Operator N 2 has eigenvalues N(N + 1) with (half)integer values of N. 

\Dm) = \ N = K -1/2, N z = m- 1/2} m^-K-1 

Thus, degenerate states with m ^ —K — 1 transform under the 2fT-dimensional representation of this 
su(2), while nondegenerate states are singlets. 

Note that by selection rules taking any other vector, e.g. a linear combination of K and S, instead 
of K in the definition of L, equation (15), will produce the same results and lead to the same su(2) 
algebra. 

As was mentioned above, to connect the state \D, —K — 1) to the rest of the degenerate subspace, 
one needs an operator that can change m by at least two. Let us consider, for example, operators 

M ± = PSlP (20) 

Using operators M± and operators iVj, one can construct a new su(2) algebra with a representation 
that incorporates the state \D, —K — 1). We define 



J + = ( J_)t = J 1 ^- 1 N + - + 3/2 PSl P_ K _, 

(21) 

Jz = G Z P + P-K-l 
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Here 3 P-k-i denotes the projector onto the state \D, —K — 1) 

P-k-i\D, -K-l) = \D, -K - 1) P- K -i\m ^ -K - 1) = 
The operator P-k-i can be written in terms of G z as follows 

= TO JL (m - 6J 

One can check by a direct computation using equations (21), (19), (17), and (14) that 

1. Operators J± and J z commute with Hamiltonian (1) at x — 1 and form an su(2) algebra. 

2. All degenerated states are uniquely labelled by eigenvalues of J z and J 

\Dm) = \J = K, J z = m) - K + 1 < m < K 

and 

\D, -K - 1) = | J = K, J z = -K) 

Thus, J z runs from — K to K with no gap and the degenerate subspace transforms under a 
(2K+l)-dimensional representation of the su(2) defined by equations (21). 

3. All nondegenerate states are singlets, i.e. they have J = J z = 

Note that operators M± are linear combinations of components of a conserved rank 2 tensor 
Mij = PSiSjP. By the same argument as for vectors, inserting any other rank 2 tensor instead of 
operator Y in equation (10) will result in equivalent operators. Higher rank tensors are not needed. 
Indeed, operators Jj already connect all degenerate states. Therefore, components of tensors of a rank 
greater than two can be written in terms of products of operators Jj in the same way as components 
of any tensor that commutes with Hamiltonian (1) at x — can be expressed through G x , G y , and 



4 Large K Limit 

To gain some insight into the conservation laws derived in the previous section let us consider the 
limit K ^> 1. In this limit the meaning of symmetry turns out to be transparent. Our intention in 
this section is not to present a rigorous analysis, but rather to develop an intuition about symmetries 
responsible for degeneracy. 

If K is a classical vector, we can write the Hamiltonian at x = 1 as 

H d = \(K + l/2)z + K\S n (22) 

3 Note that P- K -i = P-k-\P = PP-k-\P and G Z P = PG Z = PG Z P. Hence, operators J± and J z denned by 
equations (21) are indeed of the general form (3). 
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where z is a unit vector along the z-axis and operator S n is the projection of S onto the axis parallel to 
the vector (iT+l/2)z+K. For each orientation of K there are three energies corresponding to S n = +1, 
— 1 and - "top", "bottom" and "middle" levels . Middle levels have zero energy, independent on the 
direction of K. 

Now let us take into account quantization of K in the limit K 1. We assume nevertheless that 
S n still can be interpreted as a projection of S . 

The direction of K is specified by K z m m. States \S n = 0,G Z = m) are degenerate with respect 
to m. Evidently, symmetries responsible for this degeneracy "rotate K" for states with S n = 0, while 
keeping it unchanged for states with S n ^ 0. Generators for such rotations are 

Kv = (! " S 2 n )K x , y (l - Si) (23) 

Operators commute with the Hamiltonian (22) which follows from the spin-1 identity 

S n (l - Si) = (24) 

In deriving (23) we made use of a projection operator 

P d = l- Si P d \Dm) = \Dm) P d \Nm) = (25) 

where \Dm) and \Nm) denote degenerate and non-degenerate states respectively. Operators L^ y 
together with Lf — (1 — Sl)K z (l — Si) form an su(2) algebra 

[L?,Lf]=ie ijk L« + 0(l/K) (26) 

The degenerate subspace transforms under a 2K + 1 dimensional representation of this su(2), while 
nondegenerate states are singlets. Neglecting terms of the order of l/K, we can write Hamiltonian 
(22) as 

H d = B(G z )S n (27) 

Note that the operator B(G Z ) defined by equation (12) can be interpreted as an "effective magnetic 
field operator", which explains notation (12). Combining equation (27) with spin-1 identity (24), we 
obtain 

Si = , d P d = l- d (28) 
n B(G Z ) 2 B(G Z V 



Z/ 

T2 ' 



B(G Z Y - H$\ = (29) 

Now let us compare these results to the "quantum" case. First, one can check that operators 
Lf , H d , and P d defined above are indeed the limits of the corresponding operators introduced in the 
previous section. Namely, 

H^H d P^P d Ji^Ni^L^Lf asK^oo (30) 



The large K limit can be treated more accurately by introducing a preferred direction for K and representing K in 
terms of Holstein Primakoff bosons. 
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Therefore, we can think of the operator 




are the limit of equation (11). Equation (30) shows that the difference between TV, and Jj disappears 
in the limit of large K. Thus, the gap in the values of m has no analogs in the leading order of this 
limit. Note also that in the leading order in 1/K the middle levels of Hamiltonian (1) are degenerate 
at any x and for all integer values of S. An analysis of the subleading corrections to the large K limit 
may provide a simple explanation of why the degeneracies occur only for S — 1. 



In this section we deviate somewhat from the main topic of the paper and discuss generalized sym- 
metries of Hamiltonian (1). Generalized symmetries are anticommutation relations that explain, for 
example, the following relation between energies of top and bottom levels (see Section 2): 



We also show that the analogy to the large K limit can be pushed even further if we replace the 
Hamiltonian (1) at x — 1 with an effective Hamiltonian by keeping only the blocks with |m| < K — 1. 
First, we note that there exists a basis where the block of the Hamiltonian (5) at x — 1 reads 



Indeed, matrix (32) has the same eigenvalue equation as matrix (5) at x — 1. Therefore, matrices (32) 
and (5) represent the same operator in two different basises. From equation (31) we find that 1\ for 
I ml < K in this basis is 



5 Generalized Symmetries 



Exm + Ebtu — —3/2 



ml < K - 1 



-1 B(m) 
B(m) -1/2 
-1/2 



(32) 



1 
1 




(33) 



Comparing matrix (32) and equation (27), and using (33), we conclude that it is natural to define /, 



n 



as 



1 

1 




(34) 
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Thus, we can write (cf. (27)) 

H eff = B(G z )h - I/4/3 - 1/4/7 " 1/2 



(35) 



where h, I 2 , and / 3 are the first three Gell-Mann matrices of su(3). 



' 


1 


" 




" 


—% 


" 




' 1 





" 


1 








h = 


i 








h = 





-1 





































Matrices I\, I 2 , and I3 generate an su(2) subgroup of su(3) known as the isospin subgroup. We note 
that 

H eff i 2 + i 2 H eff = ~i 2 (36) 

Anticommutation relations of this type are sometimes called generalized symmetries. Equation (36) 
shows that if \m) is an eigenstate of the Hamiltonian at x — 1 with \m\ < K and the energy E, I 2 \m) 
is either also an eigenstate but with the energy —3/2 — E or I 2 \m) = 0. Examining the eigenstates of 
matrix (32) one can verify that I 2 \Nm) 7^ and I 2 \Dm) = 0. Therefore, we get 

E Tm + E Brn = ~ i 2 \Tm) = \Bm) I 2 \Bm) = \Tm) I 2 \Dm) = (37) 

Note also that since for each \m\ < K the operator K ■ S has three eigenvalues —K — 1, —1, and K 
E Tm + E Drn + E Bm = Tr(H) = (K + l/2)Tr(S z ) + Tr(K • S) = -2 (38) 



6 Hydrogen Atom 

Here we derive the nontrivial conservation laws for the Hydrogen atom from the general expression (3). 
The Hamiltonian in atomic units (see e.g. [1]) is 

«=f-i (39) 

In addition to the angular momentum L = r x p, Hamiltonian (39) also conserves the Runge-Lenz 
vector 

A = ^-i(pxL-Lxp) (40) 

Integrals of motion (40) explain the degeneracy in angular momentum. Components of vectors L and 
A can be combined to yield generators of the so(4) symmetry [2] of Hamiltonian (39). 

Let us show that conservation laws (40) can be derived directly from equation (3) using only the 
knowledge of exact eigenstates of Hamiltonian (39). The conservation of angular momentum is not 
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specific to 1/r potential. We therefore assume that we know that the angular momentum is conserved, 
but do not know any other conservation laws of Hamiltonian (39). 
According to (3) the conservation laws are of the form 

X = Y / PnYP n (41) 

n 

where P n projects out the eigenspace of (39) with eigenvalue E n = — l/2n 2 . We take a straightforward 
approach and try the simplest choices for Y. 

First, we look for a conserved scalar and try Y\ = f 2 and Y 2 — P 2 - This choice as well as any 
other powers of |p| and |r| result only in trivial commuting operators - combinations of H and L 2 . 
For example, 

Next, we search for a conserved vector and try Y — f . To evaluate 

B = ^P n rP n (42) 

n 

explicitly it is advantageous to use parabolic coordinates. The eigenstates of discrete spectrum are 
specified by three integers \ni,ri2,m), where m is the z-projection of the angular momentum and the 
principle quantum number is n = n\ + n 2 + \m\ + 1 (see e.g. [1]). Evaluating the matrix elements of 
B z in this basis, we find that only diagonal matrix elements are nonzero and 

3n 

B z \ni,n 2 ,m) = ~^( n i ~ n 2 )\n 1 ,n 2 ,m) (43) 

To express B z through components of r and p we have to identify operators that have eigenvalues 
n, ni, and n 2 . All we need to do is to follow backwards the derivation of eigenstates in parabolic 
coordinates. Substituting these operators instead of n, ni, and n 2 into equation (43), we find 

a sX 

B 7 = 



4H 

Since the choice of z-axis is arbitrary, we conclude that 

A 

B = 



AH 



with A given by equation (40). 

Thus, the substitution of Y — f into equation (3) produces the Runge-Lenz vector - the nontrivial 
conservation law responsible for the "accidental" degeneracy in Hydrogen atom. One can check that 
the alternative choice Y = p yields the same result. 
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7 Conclusion 



We have derived a complete set of symmetry generators of Hamiltonian (1) at x — 1 and established 
that the symmetry group of the Hamiltonian is SU(2). The degenerate subspace transforms under a 
(2K + l)-dimensional representation of this SU(2), nondegenerate states being singlets. For 
the degenerate states correspond to S n = where n is a unit vector along z + ~K/K and the symmetry 
operators have a meaning of rotations of K for these states. We outlined a general approach for finding 
symmetries that can be used when the eigenstates and eigenvalues of the Hamiltonian are known. The 
formal starting point in this approach is the general expression for conservation laws (3). We have seen 
that for both Hamiltonian (1) and the Hydrogen atom simplest choices for Y produce all nontrivial 
conservation laws. 

We do not have a satisfactory explanation of why multiple degeneracies at nonzero value of x occur 
only for S — 1 and not for other integer values of S. If this question is to be answered on the symmetry 
grounds, further insight into the nature and simple manifestations of symmetries derived in this paper 
is needed. In particular, the analysis of subleading corrections to large K limit may be an interesting 
avenue to pursue. 



8 Appendix: Commuting Operators 

Here we prove equation (3) for Hamiltonian (1). The proof can be readily generalized to any Hermitian 
operator whose eigenstates span the Hilbert space. 

In the energy representation Hamiltonian (1) can be written as 



H 



■\8 



V 



(A.l) 



where 8 is a (2K + 1) x (2K + 1) identity matrix and T> is a diagonal matrix with no two diagonal 
elements equal to each other or to —1/2. The block — 1/28 corresponds to the degenerate subspace, 
T> to the remaining levels. Since any matrix commutes with an identity matrix and only diagonal 
matrices commute with T>, an operator X commutes with the Hamiltonian if and only if it is of the 
form 

'AO] \ A 0] [00" 

(A.2) 

where A is an arbitrary (2K + 1) x (2K + 1) matrix and V is a diagonal matrix of an appropriate 
size. Using the projection operator P, we can rewrite (A.2) as 





' A 







' A 


" 




" 





X = 












+ 








_ 


V _ 




_ 


_ 




_ 


V _ 



X = PYP + A 



(A.3) 
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where Y is an arbitrary operator and A represents the second operator on the left hand side of equation 
(A. 2). The action of the operator A on the degenerate subspace is trivial. Namely, 

k\Dm) = 

Thus, with no loss of generality the search of commuting operators can be restricted to operators of 
the form 

X = PYP (A.4) 
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